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It is shown that the coefficients in the Laguerre expansions for successive generations 
in a renewal process are related by an algebraic convolution. Thus the calculations are 
easily mechanized for computation. 



In a paper published in 1940, Brown [1]^ studied 
the solution of the renewal equation under the 
assumption that the functions aj)p(»arino; in it could 
be represented by a Pearson type III function. The 
natural generalization of this type of solution is to 
assume that the functions involved in the renewal 
equation can be expanded in terms of a series of 
generalized Laguerre functions. Recently Zelen 
and Dannemiller [2] have presented a summary of 
results concerning the expansion of reliability func- 
tions in terms of generalized Laguerre functions, 
but did not discuss these expansions in the context 
of their application to renewal equations. In view 
of the increasing number of applications of renewal 
theory and the theory of semi-Markov processes 
[3-6] to stochastic models in relial)ility, biology, 
and the theory of traffic it is of some interest to re- 
turn once more to the subject of approximate 
solutions to the renewal equation. In this note 
we shall outline some of the results relating solutions 
of the renewal equation and expansion into series 
of Laguerre functions, without going too deeply into 
matters of uniqueness and convergence of the re- 
sulting series. 

We shall consider the simple renewal equation: 

<p{t)=f(t)+£<piT)!l(t-r)dT (1) 

and assume in what follows the existence of a unicjue 
solution to this equation. For most applications 
the f unctions /(^) and g(t) are assumed to be nonnega- 
tive over the interval 0<t<^cx> and the definite 
integrals 



/» 00 /» 

.fit)dt, 
Jo Jo 



g{t)dt 



are assumed to be bounded. A formal solution to 
eq (1) can be given in terms of tlie Laplace trans- 
forms o[J(t) and (i(t)j and provided that /(f) and g{t) 
are measurable, nonnegative, and l)Ounded in every 
interval, it can be shown that a unique nonnegative 



1 Figures in brackets indicate the literature references at the end of this paper. 



solution of the renewal equation indeed exists. 
Furthermore the asymptotic properties of (p{t) can 
be obtained from the properties of the individual 
Laplace transforms [7]. 

If we define the following Laplace transforms 



Jo 



"^{t)dt 



f*{s)= r e-^[fit)dt 
Jo 

then the transformed solution to eq (1) can be 
written 



<P*(s) = 



!-£?*(*•) 



(2) 



This solution has been used by Feller [7] to discuss 
in a rigorous manner an expansion first proposed by 
Lotka. The solution in question can be written 



^(t) = ^A,e^^' 



(3) 



where it is assumed that there is a deniunerable 
sequence of distinct roots, Sq, Sj, S2, ... to the 
equation 

g*is) = l. (4) 

Feller proves the following theorem: 

In order that the solution to eq (1) be expressible in 
the form of eq (3), where the series converges absolutely 
for t>0 it is necessary and sufficient that the Laplace 
transforrn 0*(s) admit the expansion 



1 — a*( 



g'^isY 



k 



' — St 



and that 2|Ak| converges absolutely. 
Ak are determined by 

A,= -r is,) I g*{s,). 



(5) 

The coefficients 



(6) 
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It is easy to generalize this theorem to cover the 
possibility of having multiple roots. One can see 
that the conversion of the renewal equation into 
an algebraic equation in the transform domain 
leads to a simple representation of the solution. 
We shall see that this simplification is retained when 
j{t) and g{t) are expanded in terms of generalized 
Laguerre functions. Hadwiger [8] has shown that 
if j{t) and gif) are represented by Taylor series 
in a finite interval and vanish outside of the interval 
then the renewal function can also be expanded as a 
Taylor series and the coefficients obey a discrete 
convolution relation. This result is very similar 
to the results to be presented here. 

To begin our exposition, let us recast eq (1) into 
a slightly different form by defining the successive 
generation functions [ipk{t)} by 

<Pn+i(t)-= ( <Pnir)git-T)dT,n=Q, 1,2, ... , (7) 
Jo 

The solution, (p{t), is then. 



^(0=S^«(0- 



(8) 



Clearly the (pn.(t) are also the successive terms in a 
solution to eq (1) by means of successive approxi- 
mations. If we define the Laplace transforms 

<pt{s)=^L{^,,(t)} (9) 

then the Laplace transform of eq (7) is 

<pUs)=r{s), ^,^,{s) = ^t{s)g'^{s) (10) 



or 



<p%^)=j*{-^w{'^)r 



(11) 



We shall essentially use the forms of the renewal 
equation given in this paragraph. 

The generalized Laguerre functions i^(x) will be 
defined by 



L7(t)- 






e't- 



d" 



dt 



„(e-'r+"') 






i-ty 
i'- 



(12) 



for n an integer. These functions have the 
normalization property 



j; 



e-H^L^{t)LT{t)dt = ^^'^^J^'^^ 8rs 



rl 



and their Laplace transforms take the form 



/: 



e-'"t'"L'S{t)dt-- 



(13) 



(14) 



We shall assume that the functions /(^) and g{t) can 
be expanded into series of the generic form 



vity 



o-\t 



1 (m-hlj 7i=o 



(15) 



where X is a parameter with dimensions of t~^. It 
is put in because in most applications the functions 
jp{t) are densities. The best results on the pointwise 
convergence of this series are contained in a paper 
by Uspensky [9]. He proved that the following 
conditions are sufficient for pointwise convergence 
(except at discontinuities, where the sum converges 



io^{jp{t+)+v(t- 



-)): 



/» 00 

1. r'e~^p'^{t)dt exists for a certain r. 

2. For — l<m<— 2 rer'f{t)dt exists for a 
certain /3 while for m>— ^only 



i 



t^ ''\p{t)\dt 



exists. 

3. The function p{t) is of bounded variation and 
absolutely integrable in any finite interval. In the 
cases of present interest most functions p{t) will 
satisfy these three criteria. In particular if the 
function /(i^) and g{t) are bounded over (0,oo) then 
conditions 1 and 2 will always be satisfied. Specific 
applications of this expansion have been made in 
references [11] and [12]. 

If the expansion of eq (15) is assumed to exist then 
the coefficients an are given by 



"" /m+n\ 



^ f L^(U)p(t)dt 



1 f.(rn+n\(-iy^: 

+n\pi\n-j/ j\ 



/m+n\ 
\ n ) 



where the jij are the moments 



Mj 



(16) 



Jo 



t^,-= 



.'p{t)dt. 



(17) 



We will assume that all of the functions jp{t) that 
are to be treated are such that moments of all orders 
exist. If ;p*(s) denotes the Laplace transform of 
p{t) then an application of eq (14) shows that 



p*{s)- 



vm+l 



'm+n 






(18) 



Thus we see that in the s domain the expansion of 
^*(s) is essentially a power series development in the 
variable s/(s+X). 

Now let us assume that j{t) and g{t) have the 
expansions 
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.m- 



r(m+i) t^» 



J^jJJ^iU) 



r(r+l) ^0 



S 9nLl,{\t) 



(19) 



and calculate the coefRcieiits in the expansion of 
ipj{t). Tlie simplest procedure is to use Jjaplace 
transforiiis. In this way we find 






(20) 



where the c[^^ are constants wliicli satisfy the 
recurrence relations 



i=0 \ J / 



(21) 



That is to say, the r's nuiy he found by successive 

convolutions with the sequence < ( • ) Oj r' Thus 

the couvolutiou integrals are converted into more 
convenient alg-ebraic convolutions. 

If the inverse transform of eq (11) is taken, the 
expansion of (pf,(t) ui the time donudn becomes 



^k(t)-- 



"" ( 



1 



m+k{r-\-l) + j\ 
m+k(r+l) ) 



\m+A;(?- + l) + l/w+A:(r+l) 



:m+k{r+\))\ 



Lj+k(r+l)^^^y (22) 



The iiiost useful case in practice is that for which 
m = r=0, i.e., an expansion in terms of l^aguerre 
functions of order zero. For this special case one has 



<P,(t)-- 



i=0 



(/.-) 



1 



\fc+lfk 



14 (\t). (23) 



This expression can be further reduced through the 
identity 



Lj(x)=^nr'{x) 



r = 



(24) 



lA"i^)=i2i(' ^J^l\ ^)lAx) m>l. (25) 

On substituting this expression into cq (22) and 
collecting tornis we find 



<P,it) = 



X*+'<* 



k\ 



///=0 j=in 



\ m-\ ) 



ey) 



k>\. (26) 



When g(j) is such tliat only a few of tlie coefficients 
{}m are nonzero, then the pi'ocechu'e of taking alge- 
braic convolutions is quite simple. For example, if 
we drop the restriction that /(O and (j(1[) be positive, 
and assume thaty(^) and (j{j) are given by 



/(0=a6-^^L„(X0 



(27) 



then a direct calculation leads to (he rc^siih for 



Cs'+X) 



(28) 



which implies 



^k(t) =a¥ V- ^'^ "'' 7 X 



{n+k7n+k\ 
\ k ) 



k\ 



■^n + km\A^) 






n+n+km-j- 1 \ j^ .(^^) _ (29) 



The simplest way to find (p{t) is to sum over k in the 
transform. This leads to a result analogous to that 
of Lotka in eq (3). U j\t) and (j{t) are given by 



fif) = !lit) = e-"(l+aL„{U)) 
then it is easily verified l)y induction that 



Cj*)=(^f)-f (J) CUirin+(^) CiH 



(0) I 
J-2n\ 



wliere 
hence 



(/f =cf j>0, r/f=0 ./<0 



cf>=6,,o+fl5;.« 



(30) 



(31) 



(32) 



(33) 



we see that cj^^=0 unless j is an integral multiple of 
n. A further calculation starting from eq (31) 
shows that 



.'«=(^+l)a' 



r=0, 1,2, . . .^+1. (34) 



With this explicit formula, the infiinte sum in eq (23) 
is reduced to a finite series 






(k+n)r+k-l k+l 



fr7i-\-k—m—l\ 

(^m+k\ \ r J 



^m+k\ 



(35) 



k y 
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There are several extensions of this method that 
are possible. It is not too difficult to generalize the 
solutions here to finite sets of semi-Markov equations 
(i.e., the matrix analog of eq (1)). It has also been 
pointed out by Harold Kaplan that discrete renewal 
equations can also be handled by the same methods 
provided that we start with Gottlieb polynomials 
[10]) rather than Laguerre functions for the relevant 
expansions. The obvious advantage of the present 
method is the fact that one can reduce the integral 
equation in an exact way to recurrence relations 
which may then be easily handled by a computer. 
The method is probably not too good if there are 
many nonzero coefficients or if only the moments of 
fit) and git) are known experimentally and the co- 
efficients must be calculated from them. Prelimi- 
nary experiments show that the results are very 
sensitive to errors in expressions for higher moments. 
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